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ABSTRACT 

We present the second extensive study of the coronal line variability in an active 
galaxy. Our data set for the well-studied Seyfert galaxy NGC 5548 consists of five 
epochs of quasi-simultaneous optical and near-infrared spectroscopy spanning a period 
of about five years and three epochs of X-ray spectroscopy overlapping in time with it. 
Whereas the broad emission lines and hot dust emission varied only moderately, the 
coronal lines varied strongly. However, the observed high variability is mainly due to 
a flux decrease. Using the optical [Fe vii] and X-ray O vii emission lines we estimate 
that the coronal line gas has a relatively low density of Ue ^ 10^ cm“^ and a relatively 
high ionisation parameter of logC/ ~ 1. The resultant distance of the coronal line gas 
from the ionising source of about eight light years places this region well beyond the 
hot inner face of the dusty torus. These results imply that the coronal line region is 
an independent entity. We find again support for the X-ray heated wind scenario of 
Pier & Voit; the increased ionising radiation that heats the dusty torus also increases 
the cooling efficiency of the coronal line gas, most likely due to a stronger adiabatic 
expansion. The much stronger coronal line variability of NGC 5548 relative to that 
of NGC 4f5f can also be explained within this picture. NGC 5548 has much stronger 
coronal lines relative to the low ionisation lines than NGC 4f 5f indicating a stronger 
wind, in which case a stronger adiabatic expansion of the gas and so fading of the line 
emission is expected. 

Key words: galaxies: Seyfert - infrared: galaxies - X-rays: galaxies - quasars: emis¬ 
sion lines - quasars: individual: NGC 5548 


1 INTRODUCTION 

In addition to the broad and narrow emission lines, the spec¬ 
tra of active galactic nuclei (AGN) display high-ionisation 
emission lines, the so-called coronal lines, which require en¬ 
ergies ^ 100 eV to be excited. The coronal line region is 
believed to lie at distances from the central ionising source 
intermediate between those of the broad (BELR) and nar¬ 
row emission line region (NELR) and to possibly coincide 
with the hot inner face of the c ircumnuclear, obs curing 
dusty torus (as first suggested bv IPier fc VoitiriOOSl l. This 
assumption is based mainly on the following: (i) the coro- 
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nal lines have 3 — 4 orders of magnitude higher critical 
densities for collisional deexcitation than the low-ionisation 
narrow emission lines and are believed to be emitted at 
these densities; (ii) the coronal line profiles often but not 
always have full widths at half maxima (FWHM) inter¬ 
mediate between those of the broad and n arrow emission 
lines (FWHM^ 500 — 1500 km s~^; e.g.. IPen ston et ^ 
]j84_ ApT3enzeller_fcOestreicherl 19881 : lErkens et al.l llQO*^ 
Rodriguez-Ardila et al. 2002I . 201lll : and (iii) the emission 
from this region is often extended but much less so than 
that from t he low-ionisation NELR (on scales of ~ 80 — 
150 pc; e.g..|Prie to et al.l l2005l: iM iiller Sai ichez et ^l2006l : 


Muller^Sanchez_etal. 20111 '' Mazzalav et al.l~13l b Recently. 

Rose et al.l ( 2015h found evidence in their study of seven 
AGN with strong coronal lines selected from the Sloan Dig¬ 
ital Sky Survey (SDSS) that the coronal line region has 
higher densities than the low-ionisation NELR and lies at 
distances from the ionising source similar to those estimated 
for the hot dust sublimation radius. 
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In iLandt et alJ (l2015l ~l , we have presented the first ex¬ 
tensive study of the coronal line variability in an AGN and 
have started to question this picture. Our data set for the 
nearby, well-known source NGC 4151 included six epochs of 
quasi-simultaneous optical and near-IR spectroscopy span¬ 
ning a period of ~ 8 years and five epochs of X-ray spec¬ 
troscopy overlapping in time with it, with the observations 
in each wavelength performed with the same telescope and 
set-up. For this source we found that the coronal lines varied 
only weakly, if at all, and estimated a relatively low density 
(rie ~ 10® cm“®) and high ionisation parameter (log 17 ~ 1) 
for the coronal line gas. The resultant distance of the coro¬ 
nal line region from the ionising continuum was ~ 2 light 
years, which is well beyond the hot inner face of the ob - 
scuring torus of (~ 2 light months; iKoshida et ahlEoidll . 
Based on the high ionisation parameter, we proposed that 
the coronal line region is an independent entity rather than 
part of a continuous gas distribution connecting the BELR 
and low-ionis ation NELR , possi bly an X-ray heated wind as 
suggested by IPier fc VoitI dlQQSh . 

The most stringent constraints on the properties of 
the coronal line emitting region could come from variabil¬ 
ity studies, in particular if the variability of several coronal 
lines can be compared with each other and with that of other 
AGN components such as the BELR, h ot dust emiss i on an d 
X-ray luminosity, as we have done in iLandt et al.l ll2015h . 
However, mainly due to the weakness of these emission lines 
and also lack of d ata, very f ew stu dies of this kind have been 
attempted so far. IVeilleu:^ lll988ll made the only systematic 
study of the coronal line variability. In his sample of ~ 20 
AGN he found firm evidence that both the [Fe VIl] A6087 
and [Fe x] A6375 emission lines varied (during a period of 
a few years) for only one source, namely, NGC 5548, which 
is the subject of this paper, and tentative evidence for an¬ 
other seven sources (including NGC 4151). Then, within a 
general optical variabilit y campaign on the sourc e Mrk 110 
lasting for half a year, iKollatschnv et ahl (|200ll l reported 
strong [Fe x] variations. More recently, follow-up optical 
spectroscopy of a handful of objects with unusually promi¬ 
nent coronal lines selected from the SDSS showed that in 
half of them the coronal lines strongly faded (by factors of 
~ 2 — 10), making these sources candidates for stellar tid al 
disruption events llKomossa et al.ll2009l : lYang et al.ll2013fl . 

Here we present the second extensive study of the coro¬ 
nal line variability in an AGN. Our data set for the well- 
studied source NGC 5548 {z = 0.017) consists of five epochs 
of quasi-simultaneous optical and near-infrared spectroscopy 
spanning a period of ~ 5 years and three epochs of X-ray 
spectroscopy overlapping in time with it. The paper is or¬ 
ganised as follows. In Section 2, we present the data and 
measurements. In Section 3, we discuss the observed vari¬ 
ability behaviour of the near-IR, optical and X-ray coronal 
lines, for which we seek an interpretation in the context of 
the location of and excitation mechanism for the coronal 
line emission region in Section 4. Finally, in Section 5, we 
summarise our main results and present our conclusions. 
Throughout this paper we have assumed cosmological pa¬ 
rameters Ho = 70 km s“® Mpc“®, Hm = 0.3, and Ha = 0.7. 


2 THE DATA AND MEASUREMENTS 

2.1 Near-IR and optical spectroscopy 

We have five epochs of quasi-simultaneous (within less than 
two months) near-IR and optical spectroscopy for NGG 5548 
(see Tables [T] and [2]) . All data are of relatively high signal- 
to-noise ratio (continuum S/N > 50 — 100). The near- 
IR spectroscopy wa s obtained with the SpeX spectrograph 
dRavner et al.[|2003l i at the NASA Infrared Telescope Fa¬ 
cility (IRTF), a 3 m telescope on Mauna Kea, Hawai’i, in 
the short cross-dispersed mode (SXD, 0.8 — 2.4 /rm). All 
data were obtained through a slit of 0.8 x 15” giving an 
average spectral resolution of full width at half maximum 
(FWHM) ~ 400 km s“^. The four epochs spanning the years 
2004 — 20 07 are our own dat a an d were discus s ed an d pre¬ 
sented in ILandt et al.1 ll2008h and ILandt et al.1 |201_l|). The 
near-I R spectrum from 2002 was discussed by iRiffel et al.1 
1I2OO6II . The o ptical spectra were o btained with the FAST 
spectrograph dFabricant et al.lll998ll at the Tillinghast 1.5 m 
telescope on Mt. Hopkins, Arizona, using the 300 1/mm grat¬ 
ing and a 3" long-slit. This set-up resulted in a wavelength 
coverage of ~ 3720 — 7515 A and an average spectral res¬ 
olution of FWHM ~ 330 km s“®. The slit was rotated to 
the parallactic angle for all observations but the May 2004 
epoch. However, the latter spectrum was observed at a very 
low air mass (sec 2 1.01). The optica l da ta were dis¬ 

cussed and presented in ILandt et al.l (I2OO8II and lLandt et al.1 
ll201lll . with the exception of the May 2004 spectrum, which 
was retrieved from the FAST archive. 

We have measured the fluxes of the strongest near- 
IR and optical coronal lines and es t imate d their la uncer¬ 
tainties as detailed in ILandt et al.1 (l2015l L In the near-IR, 
we have measured two sulfur lines and two silicon lines, 
namely, [S VIIl] A9911, [S ix] 1.252/rm, [Si Vl] 1.965pm 
and [Si x] 1.430pm (see Table (TJ. In the optical, we have 
measured four iron emission lines, namely, [Fe VIl] A3759, 
[Fe VIl] A5159, [Fe Vii] A5721 and [Fe Vii] A6087 (see Ta¬ 
ble [21l. The Icr uncertainties are ^ 1 — 5% for the strongest 
lines and ~ 2 — 20% for the weakest ones. Since the spectra 
were obtained in non-photometric sky conditions, we study 
in the following the temporal changes of the coronal lines 
in relative rather than absolute flux. In particular, we scale 
the coronal line emission to that of a strong, forbidden low- 
ionisation emission line that is unblended and observed in 
the same spectrum. The emission region that produces the 
low-ionisation narrow lines is believed to be located a t large 
enough distances from the central io nising source (e.g. lPoera 
ll988lJThdhunter fc Tsvetanovlll989h and to have a relatively 
low density for its flux to remain constant on timescales of 
decades. We have scaled the near-IR and optical coronal 
lines to the [S ill] A9531 and [O ill] A5007 emission lines, 
respectively. 

In order to further constrain the significance of the ob¬ 
served coronal line variability, we have considered the two 
extreme cases of where no variability and the highest vari¬ 
ability are expected. For both scaling lines we observe also 
the other emission line that is emitted from the same up¬ 
per level, namely, [S III] A9069 and [O ill] A4959. Their ob¬ 
served ratios, which should be close to the theoretical values 
of [S III] A9531/A906 9=2.58 and [O ill] A5007/A4959=2.92 
llKramida et al.lr2013h . are not expected to vary and so their 
observed variability sets a lower threshold for the signlfi- 
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Table 1. Near-IR emission line fluxes and ratios 


Observation 

[S III] A9531 

[S III] A9069 

[S III]/ 

[S VIII] A9911 

[S III]/ 

[S IX] 1.252/rm 

[S III]/ 

Date 

(erg/s/cm^) 

(erg/s/cm^) 

[S III] 

(erg/s/cm^) 

[S VIII] 

(erg/s/cm^) 

[S IX] 

2002 Apr 23 

(2.79±0.06)e-14 

(1.24±0.04)e-14 

2.25±0.09 

(4.45±0.17)e-15 

6.3±0.3 

(2.88±0.17)e-15 

9.7±0.6 

2004 May 23 

(4.50±0.11)e-14 

J.76±0.11)e-14 

2.56±0.17 

(4.02±0.76)e-15 

11.2±2.1 

(2.60±0.99)e-15 

17.3±6.6 

2006 Jan 9 

(3.86±0.22)e-14 

J.79±0.11)e-14 

2.16±0.18 

(5.22±0.73)e-15 

7.4±1.1 

(4.44±0.93)e-15 

8.7±1.9 

2006 Jun 12 

(6.49±0.15)e-14 

(2.65±0.19)e-14 

2.45±0.19 

(4.89±0.55)e-15 

13.3±1.5 

(3.79±0.55)e-15 

17.1±2.5 

2007 Jan 24 

(5.36±0.10)e-14 

(1.96±0.07)e-14 

2.74±0.11 

(4.23±0.53)e-15 

12.7±1.6 

(2.73±0.29)e-15 

19.6±2.1 

Observation 

[S III] A9531 

[Si VI] 1.965/rm 

[S III]/ 

[Si X] 1.430/j-m 

[S III]/ 

Pa/3 

[S III]/ 

Date 

(erg/s/cm^) 

(erg/s/cm^) 

[Si VI] 

(erg/s/cm^) 

[Six] 

(erg/s/cm^) 

Pa/3 

2002 Apr 23 

(2.79±0.06)e-14 

(1.42±0.02)e-14 

2.0±0.1 

(6.57±0.21)e-15 

4.3±0.2 

(6.90±0.10)e-14 

0.404±0.011 

2004 May 23 

(4.50±0.11)e-14 

(1.30±0.18)e-14 

3.5±0.5 

G.19±0.12)e-14 

3.8±0.4 

J.ll±0.04)e-13 

0.405±0.017 

2006 Jan 9 

(3.86±0.22)e-14 

J.32±0.05)e-14 

2.9±0.2 

(3.38±0.55)e-15 

11.4±2.0 

J.26±0.03)e-13 

0.306±0.020 

2006 Jun 12 

(6.49±0.15)e-14 

(2.27±0.07)e-14 

2.9±0.1 

(5.86±0.45)e-15 

11.1±0.9 

J.68±0.05)e-13 

0.386±0.015 

2007 Jan 24 

(5.36±0.10)e-14 

J.01±0.04)e-14 

5.3±0.2 

(3.61±0.23)e-15 

14.9±1.0 

J.37±0.02)e-13 

0.391±0.010 


Table 2. Optical emission line fluxes and ratios 

Observation 

Date 

[O III] A5007 
(erg/s/cm^) 

[O III] A4959 
(erg/s/cm^) 

[O III]/ 

[O III] 

[Fe VII] A3759 
(erg/s/cm^) 

[O III]/ 
[Fe VII] 

[Fe VII] A5159 
(erg/s/cm^) 

[O III]/ 

[Fe VII] 

2002 Mar 6 

(4.72±0.06)e-13 

(1.51±0.06)e-13 

3.13±0.13 

(1.94±0.24)e-14 

24.3±3.1 

(l.ll±0.20)e-14 

42.5±7.5 

2004 May 24 

(5.21±0.05)e-13 

J.90±0.07)e-13 

2.74±0.10 

J.94±0.20)e-14 

26.9±2.8 

J.29±0.14)e-14 

40.4±4.4 

2006 Jan 9 

(5.89±0.02)e-13 

(2.01±0.03)e-13 

2.93±0.05 

J.37±0.13)e-14 

43.0±4.1 

(8.31±1.01)e-15 

70.9±8.6 

2006 Jun 20 

(4.19±0.04)e-13 

(1.40±0.02)e-13 

2.99±0.06 

J.85±0.19)e-14 

22.7±2.3 

(9.03±1.00)e-15 

46.4±5.2 

2007 Feb 17 

(2.97±0.03)e-13 

(9.71±0.24)e-14 

3.06±0.08 

(8.76±1.24)e-15 

33.9±4.8 

(5.04±0.71)e-15 

58.9±8.3 

Observation 

[O III] A5007 

[Fe VII] A5721 

[O III]/ 

[Fe VII] A6087 

[O III]/ 

Ha 

[O III]/ 

Date 

(erg/s/cm^) 

(erg/s/cm^) 

[Fe VII] 

(erg/s/cm^) 

[Fe VII] 

(erg/s/cm^) 

Ha 

2002 Mar 6 

(4.72±0.06)e-13 

(1.63±0.15)e-14 

29.0±2.7 

(3.19±0.33)e-14 

14.8±1.5 

(1.77±0.02)e-12 

0.267±0.004 

2004 May 24 

(5.21±0.05)e-13 

G-49±0.15)e-14 

35.0±3.5 

(2.00±0.13)e-14 

26.1±1.8 

(2.00±0.04)e-12 

0.261±0.006 

2006 Jan 9 

(5.89±0.02)e-13 

J.02±0.06)e-14 

57.8±3.6 

J.70±0.07)e-14 

34.7±1.4 

(1.57±0.02)e-12 

0.375±0.005 

2006 Jun 20 

(4.19±0.04)e-13 

(7.30±0.59)e-15 

57.4±4.7 

J.65±0.10)e-14 

25.4±1.5 

J.30±0.02)e-12 

0.322±0.005 

2007 Feb 17 

(2.97±0.03)e-13 

(4.29±0.43)e-15 

69.2±7.0 

(8.04±0.90)e-15 

36.9±4.1 

(1.04±0.01)e-12 

0.286±0.005 


cance of the coronal line variability. Then, we have measured 
the fluxes of the two prominent broad emission lines Pa/3 
(in the near-IR) and Ha (in the optical). Since the BELR 
is expected to be the most variable of any AGN emission 
line region, the observed flux increase of these broad lines 
gives an estimate of the maximum value that can be reached 
within the current data set. These results are also listed in 
Tables [T] and [21 

In addition to the emission lines, we have measured in 
the near-IR spectra the continuum fluxes at the rest-frame 
wavelengths of 1 /xm and ~ 2.1 /rm (see Table [3]). As we 
have shown in lLandt et al.l ll201ll ~l. the former is dominated 
by the accretion disc flux, which is believed to be the main 
source of ionising radiation in AGN and so the driver of 
the observed variability, whereas the latter is emitted from 
the hot dust component of the obscuring torus. Further¬ 
more, we have derived the hot dust temperature from black- 
body fits to t he ne ar-IR spectral continuum as described in 
iLandt et al.l 1 I 2 OIII ') and list it also in Table O In Section 21 
we will compare these values to the optical coronal line ra¬ 
tios [Fe vii] A6087/A3759 and [Fe Vii] A5159/A6087, which 
are suitable indicators of the gas temperature and de nsity 
dNussbaumer et all Il982l : iKeenan k, Norringtonl Il987tl , re¬ 
spectively (listed in Table (Sj. 


2.2 X-ray spectroscopy 

As one of the brightest Seyfert 1 galaxies at X-ray fre¬ 
quencies, NGC 5548 was intensively studied in the years 
1999-2002, after which period the source became unobserv¬ 
able by XMM-Newton and had only limited observability 
with Chandra. Due to the limited observability, there are 
two shorter Chandra observations in 2005 and 2007, until 
a long observational campaign took place in the summer of 
2013 and winter of 2013/14. Between 2002 and 2007, the 
time period overlapping with our near-IR and optical spec¬ 
troscopy, there are a total of three Chandra observations 
of NGC 5548 with the low-energy transmission grating in 
combination with the high-resolution camera spectrograph 
(LETGS). The 2002 observations showed the source in an 
average flux state, while the 2005 and especially the 2007 
spectra were obtained in low-flux states. The exposure times 
were 340 ks, 141 ks and 162 ks in 2002, 2005, and 2007, re¬ 
spectively. All three X-ray observations were taken over a 
period of three days during two consecutive orbits. The low 
flux state helps in the detection of emission lines but the 
shorter exposure times lead to low S/N spectra. 

The X-ray s pectra obtained between 2002 and 2007 
were analyse d bv lsteenbrugge et al.l (120051 ') , jPetmers et al.l 
(I 2 OO 8 I . I 2 OO 9 I ) and, most recently, Ebrero et al. (2015, in 
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Table 3. Near-IR continuum fluxes and optical coronal line ratios 


Observation 

Date 

(erg/s/cm^) 

A/l,jm/ 

[S III] 

(erg/s/cm^) 

A/2.1^m/ 

[S III] 

Tdust 

(K) 

[Fe VII] A6087/ 
[Fe VII] A3759 

[Fe VII] A5159/ 
[Fe VII] A6087 

2002 Apr 23 

(7.74±0.08)e-12 

277.4±6.7 

(1.38±0.01)e-ll 

494.6±10.9 

1371 

1.644±0.268 

0.348±0.071 

2004 May 23 

(1.48±0.07)e-ll 

328.9±17.1 

(1.85±0.02)e-ll 

411.1±11.5 

1572 

1.031±0.129 

0.645±0.082 

2006 Jan 9 

(2.38±0.09)e-ll 

616.6±43.2 

(2.68±0.04)e-ll 

694.3±41.7 

1730 

1.241±0.128 

0.489±0.062 

2006 Jun 12 

(2.32±0.06)e-ll 

357.5±11.8 

(3.56±0.02)e-ll 

548.5±13.2 

1547 

0.892±0.105 

0.547±0.069 

2007 Jan 24 

(2.10±0.04)e-ll 

391.8±11.0 

(2.19±0.02)e-ll 

408.6±8.2 

1291 

0.918±0.166 

0.627±0.112 


Table 4. X-ray emission line and continuum fluxes 


Observation O VII f 0.561 keV* /o. 2 — lOkeV 

Date Ebrero et al. Detmers et al. Ebrero et al. 

(ph/s/cm^) (ph/s/cm^) (erg/s/cm^) 

2002 Jan 18 (0.83±0.07)e-04 (0.75±0.07)e-04 (1.51±0.02)e-10 

2005 Apr 15 (0.42±0.07)e-04 (0.35±0.06)e-04 {1.90±0.06)e-ll 

2007 Aug 14 (0.34±0.07)e-04 (0.27±0.06)e-04 (1.05±0.06)e-ll 

* Ebrero et al. contrary to Detmers et al. assume that the narrow emission 
lines are absorbed by the warm absorber. 


prep.). Both lPetmers et al.l (l2009l 'l and Ebrero et al. analyse 
all previously observed X-ray spectra using the latest atomic 
data. The difference between the two is that Ebrero et al. 
focus on detecting changes in the wa rm absorber with con - 
tinuum flux level and so, contrary to lPetmers et all ll2009ll . 
assume that the narrow emission lines are absorbed by the 
warm absorber. Generally, it has been assumed that the 
warm absorber does not affect the flux of the narrow emis¬ 
sion lines, however, Whewell et al. (2015, submitted) have 
recently presented evidence that at least one and potentially 
up to three warm absorber components do indeed absorb at 
least the O VII triplet and most likely all narrow X-ray emis¬ 
sion lines. In Table [d] we list the measurements and their la 
errors from both studies, Ebrero et al. and Petmers et al., 
which yield similar results. Pue to the still rather high con¬ 
tinuum level in 2005 and 2007 and lower S/N in 2005, only 
the O VII f narrow emission line was detected in all three 
spectra. 

3 THE VARIABILITY BEHAVIOUR 

3.1 The near-IR and optical coronal lines 

Over the ~ 5 year period sampled by the near-IR and optical 
spectroscopy the ionising flux of the accretion disc, which is 
assumed to be the main driver for the variability of the broad 
line, coronal line and hot dust emission regions, stayed con¬ 
stant in the first two years, increased by a factor of 2 in the 
following two years and decreased again to its original flux 
state during the last year (see Fig.[T|. We note that the flux 
increase of ~ 19% observed between the April 2002 and May 
2004 epochs is consistent with zero, once the deviations from 
the theoretical value of the [S ill] A9531/[S ill] A9069 ratios 
measured in these spectra are taken into account (~ 14.8% 
and ~ 1.1%, respectively). The variability of the ionising 
continuum is mirrored in that of the hot dust emission and 
Pa/3 broad line region (see Figs. [T] and (2] respectively); also 
their flux stayed constant in the first two years, increased 



1000 2000 3000 1000 2000 3000 

MJD (51600 I ) 

Figure 1. The variability of the continuum fluxes at rest-frame 
wavelengths of ~ 1 /rm (sampling the accretion disc), ~ 2.1 /rm 
(sampling the hot dust emission) and in the energy range 0.2 — 
10 keV (related to the accretion disc) and that of the O VII f 
X-ray coronal line. The near-IR continuum fluxes were divided 
by the (S III] A9531 line flux and normalised to the value of the 
May 2004 epoch. The X-ray measurements were normalised to 
the value of the April 2005 epoch. Note the different y-axis scale 
in the bottom two panels. We plot Icr error bars. 

in the following two years (by 70% and ~ 30%, respec¬ 
tively), and then continuously decreased to its original level 
for the next year. The flux of the Ha broad line region also 
stayed constant in the first two years, but then decreased in 
the next two years (by ~ 40%), followed by an increase to 
its original level during the last year. 

The variability response of the near-IR coronal lines 
was stronger than that of the hot dust emission and broad 
emission line region and in general very different from it. 
Whereas the flux of the ]S ix] and [Si x] lines, the lines with 
the highest ionisation potentials, stayed constant in the first 
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Figure 2. The variability of the near-IR coronal lines in the pe¬ 
riod April 2002 - January 2007. For comparison, we show also 
the [S III] A9069 narrow line and Pa/3 broad line, which are 
expected not to vary and to vary maximally, respectively. The 
[S III] A9531/]S III] A9069 flux ratios were normalised to the the¬ 
oretical value, whereas those between (S III] A9531 and the other 
lines were normalised to the value of the May 2004 epoch. We 
plot la error bars. The ionisation potentials of the coronal lines 
are given at the top right. 


Figure 3. Same as Fig. [2] for the optical coronal lines. For com¬ 
parison, we show also the [O III] A4959 narrow line and Ha broad 
line, which are expected not to vary and to vary maximally, re¬ 
spectively. 


a flux decrease by ~ 60% followed by a flux increase (by 
~ 90%) and then flux decrease (by ~ 50%). 


two years, the flux of the other two near-IR coronal lines 
decreased by a factor of ~ 2. After this initial decrease, the 
flux of the [S VIIl] line stayed constant, whereas the flux of 
the [Si Vl] line stayed constant in the following two years, 
but decreased further in the last year, resulting in a total 
decrease during the observing period of a factor of ~ 3. The 
flux of the [S ix] line decreased significantly (by a factor 
of ~ 2) only during the last year of the observing period, 
whereas the [Si x] line started to diminish already in May 
2004 and by the end of the observing period reached a flux 
decrease of a factor of ~ 4. 

Similarly to the [S VIIl] and [Si Vl] near-IR coronal lines, 
the flux of the optical coronal line [Fe VIl] A6087 decreased 
by ~ 80% in the first two years, whereas the flux of the 
other three [Fe VIl] lines stayed constant. The flux of the 
[Fe VIl] A6087 line further decreased in the following two 
years (by ~ 30%), then increased and decreased again, thus 
reaching a flux decrease of a factor of ~ 2 by the end of 
the observing period. Similarly to the Ha broad line, the 
flux of the [Fe VIl] A5159 line, after staying constant for 
the first two years, decreased in the following two years (by 
~ 80%), followed by an increase to its original level during 
the last year. Similarly to the [Si x] near-IR coronal line, 
the flux of the [Fe VIl] A5721 line started to diminish in 
May 2004 and by the end of the observing period reached a 
flux decrease of a factor of ~ 2. Finally, after May 2004, the 
[Fe VIl] A3759 line varied similarly to the [Fe VIl] A6087 line; 


3.2 The X-ray coronal lines 

In the ~ 5 year period sampled by the X-ray spectroscopy 
the unabsorbed 0.2 — 10 keV continuum flux, which is as¬ 
sumed to be produced by the central ionising source and so 
to be linked to the accretion disc flux, decreased by a factor 
of ~ 8 in the first three years, and then decreased further 
by a factor of ~ 2 in the next two years (see Fig. [T]). 

The variability response of the forbidden X-ray coronal 
line O VII f was similar to that of the near-IR and optical 
coronal lines (see Fig. (TJ. Its line flux decreased in the first 
three years by a factor of ~ 2 and stayed constant in the 
next two years. We note that this result is independent of 
the assumption of whether the X-ray narrow lines are af¬ 
fected by the warm absorber, since the ratio between the 
absorbed and un absorbed line flux is similar for all three 
observing epochs. iDetmers et al.l (l2008ll list also the flux or 
non-detection limit for the Ne IX f line, which is detected 
only in the 2002 spectrum. The trend for this much weaker 
X-ray emission line is similar to that observed for the O VII f 
line; a decrease in flux during the observing period by a fac¬ 
tor of ^ 3. Finally, it is worth noting that the flux of both 
the O VII f line and the unabsorbed X-ray continuum recov¬ 
ered by the summer o f 2013 to a level on ly slightly below 
that observed in 2002 (iKaastra et al.l[2014h . 
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4 THE ORIGIN OF THE CORONAL LINE 

EMISSION REGION 

The current understanding is that the coronal line region 
in AGN is dust-free, photoionised and located beyond the 
BELR at distances from the central ionising sonrce simi¬ 
lar to th ose of the hot inne r face of the obscnring dnsty 
torus. In iLandt et al.l (l2015fl we confirmed for the source 
NGC 4151 that the coronal line region is photoionised but 
showed that it is located well beyond the hot dusty torus, 
at a distance from the central ionising source of a few light- 
years. This explained the low variability amplitude observed 
for the coronal lines in this source. 

The situation in the source NGC 5548 seems at first 
glance to be very different. Whereas the broad emission lines 
changed only little (by ~ 30 — 40%) in the period covered 
by our data and the hot dust emission changed moderately 
(by ~ 70%), the coronal lines varied strongly (by factors of 
2 — 4). However, the observed high variability amplitude 
of the coronal lines is mainly due to a flux decrease. In fact, 
none of the coronal lines reacted with a flux increase to the 
increase in ionising flux observed between the May 2004 and 
January 2006 epochs as the Pa(l broad line and hot dust 
emission did. 

In the following, we show that the coronal line region in 
NGC 5548 is similar to that in NGC 4151, i.e. it is a separate 
entity, most likely an X-ray heated wind, that has a low 
density and a high ionisation parameter and that is located 
far away from the dusty torus fSection 14.11) . However, this 
wind appears to be stronger in NGC 5548 than in NGC 4151 
and so, most likely, it undergoes a more severe adiabatic 
expansion and cooling after being launched (Section 14.211 . 
Finally, in Section [4.31 we discuss the source NGC 5548 in 
the context of stellar tidal disruption events, which have 
been proposed as an alternative explanation for the fact that 
strong coronal line emitters often show a strong fading of 
their coronal lines over a period of several years. 


4.1 Plasma diagnostics and the similarity between 
NGC 5548 and NGC 4151 


As in iLandt et all (l2015l l. we estimate the density of 
the coronal line gas using the three optical iron lines 
[Fe VII] A3759, [Fe Vii] A5159 and [Fe Vii] A6087, 
since, as discussed by previous studies, the line ratios 
[Fe vii] A6087/A3759 and [Fe Vii] A5159/A6087 are suit- 
able indicators of temperature and density, respec tively 
llNussbaumer et ahl Il98d : iKeenan k, Norringtonl Il987l l . For 
these emission lines, we have generated a temperature ver¬ 
sus density grid using version 13.03 o f the plasma simula - 
tion code CLOUDY (last described by iFerland et al.1 l2013l l 
(see Fig.|4ll. For the photoionisation s imulations we have as - 
sumed solar abundances as given bv iGrevesse et ahl (l2010l l 
and approximated the incident radiation field with the 
me an AGN spectral energy distributio n (SEP) derived 
by Mathews fc Ferlandl (Il987l l . Recently, iMebdipour et al.1 
ll2015ll observed the hard X-ray portion of the SED and 
found that it continues to rise, in vfv, until at least ~ 
80 keV. So NGC 5548 was in a very hard state when their 
observations were done. Very hard X-rays actually have only 
modest effects on the physical conditions within a pho¬ 
toionised cloud due to their small photoelectric cross sec- 



[Fe VII] 5159/6087 

Figure 4. The observed optical coronal line ratios 
[Fe VII] A6087/A3759 versus [Fe VIl] A5159/A6087 for NGC 5548 
(filled circles) and NGC 4151 (open circles), overlaid with curves 
of constant temperature (logT = 4,4.1, 4.2, ...5 K) and constant 
number density (logrie = 3, 3.2, 3.4, ...7 cm“®) for the case of 
photoionisation equilibrium. The case of collisional ionisation 
equilibrium is shown in green. We plot la error bars. 


tion. The emission line spectrum is mainly affected by the 
FUV/XUV part of the SED, which is not dir ectly observable 
excep t in the very highest redshift objects llGollinson et al.1 
1201511 . The FUV/XUV part of the older SED was derived 
from line observations which were sensitive to that portion 
of the spectrum. The two SEDs are in qualitative agreement 
in this part of the spectrum, so our results would not have 
changed had we used the newest observations. We have as¬ 
sumed the cases of either photoionisation or collisional ion¬ 
isation equilibrium. 

The measurements for NGC 5548 (filled circles) give 
a similar answer to those for NGC 4151 (open circles); 
the coronal line gas is photoionised rather than collision- 
ally ionised and its density appears to be relatively low. 
All but one observing epochs for NGC 5548 constrain it to 
Tie ~ 10^ cm“® within ~ 3 ct, whereas the May 2004 epoch 
reaches this value within ~ 4a. We note that the estimated 
density is several orders of magnitude below the critical den¬ 
sity of [Fe vii] of ncrit ~ 3 x 10^ cm“®. We find again for 
NGC 5548, as we did previously for NGC 4151, that all mea¬ 
sured [Fe vii] ratios lie to the right of the theoretical grids. 
This indicates most likely a large ionisation parameter, since 
in such a case the contribution from fluorescence increases 
significantly the [Fe VIl] A5159/A6087 ratio (see Fig. 7 in 
ILandt et ^ 20151). 

As in Landt et al. I iliolsl) , we estimate the ionisation pa¬ 
rameter of the coronal line gas using the three X-ray lines 
from the helium-like ion of oxygen O VII f, O VII i and 
O VII r. Using the plasma simulation code CLOUDY, we have 
generated a ionisation parameter versus density grid for the 
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Figure 5. The observed X-ray coronal line ratios O VII (f-|-i)/r (G 
ratio) versus O VII f/i (R ratio) for NGC 5548 (filled symbols) 
and NGC 4151 (open circles), overlaid with curves of constant 
ionisation parameter (log C/ = —1, —0.5, 0, ...2) for two constant 
number densities (rie = 10^ cm~^ and rie = 10^^ cm“^) and a 
fixed column density of Nn = 10^^ cm“^, assuming photoion¬ 
isation equilibrium. For NGC 5548 we plot both the absorbed 
(filled circle) and unabsorbed case (filled triangle). The case of 
collisional ionisation equilibrium is shown in green. We plot Icr 
error bars. 

O VII lines (see Fig. El). As discussed by previous studies, the 
X-ray line ratios O VII (f+i)/r (the so-called G ratio) and 
O VII f/i (the so-called R ratio) trace the ionisation param¬ 
eter, which is directly related to the kinetic gas tempera- 
ture, and density, respec t ively, for a given colum n density 
JPorauet fc Dubau I 2 OOOI : iPorter fc Ferlandl l2007l l. We as¬ 
sumed again the cases of either photoionisation or collisional 
ionisation equilibrium. Neither of the studies mentioned in 
Section EJ detects all three O VII lines and, therefore, no 
density and/or ionisation parameter can be determined from 
this triplet for the period 2002-2007. However, the entire 
O VII triplet is detected bv iKaastra et all (120141 ') when the 
source was in an obscured s tate similar to what is gener¬ 
ally observed for NGC 4151. Ii5iastra et al.l (l2014l l assumed 
that the narrow emission lines are not absorbed by the 
warm absorber (similar to the study of lPetmers et al]l2009l l 
and from their data we determine a G ratio of 4.40±0.57 
and an R ratio of 3.49±0.35 (filled triangle in Fig. El). 
lAndrade- Velazquez et all ll2010h combined all Chandra spec¬ 
tra taken before 2007 and measure the O VII triplet in the 
stacked HETG and LETG spectra, also assuming that the 
lines are not absorbed by the warm absorber. From their 
data we calculate a G ratio of and d.Ojlj ® and an 

R ratio of S.Ojlj ® and 2.8lj g for the HETG and LETG 
data, respectively. Given the large errors, th ese results are 
consis tent with the XMM-Newton results of iKaastra et al.l 
(l2014h . Therefore, we will use here the latter data given the 
higher quality, although the observations are not contempo¬ 


raneous with our optical spectroscopy. If we assume that the 
O VII triplet is affected by the warm absorber, both the G 
and R ratios decrease, because the resonance and intercom- 
bination lines are absorbed whereas the forbidden line is not. 
iMehdipour et af] (l2015l ) have recently pointed out that inner 
shell lines of O VI can absorb components of the O VII lines 
used as a di agnostic here and in observations of NGC 4151 
presented in iLandt et al.l (l2015lL This hypothesis is attrac¬ 
tive because it is hard to think of any other way to account 
for the unusual ratios found in our previous paper. But the 
problem with their explanation in the case of NGC 5548 is 
that the measured O VI column density (A^ovi = 10^® cm“^; 
ISteenbrugge et 51l2005ll of the warm absorber, which po¬ 
tentially can be cospatial with the O VII line emitting re¬ 
gion, is ~ 1.5 dex smaller than that needed to produce 
the absorption assuming their atomic data. Then, using the 
line strengths from Whewell et al. (2015, submitted), the G 
and R ratios become 3.40±0.52 and 3.70±0.27, respectively 
(filled circle in Fig. El). The absorbed values constrain the 
ionisation parameter of the coronal line gas in NGC 5548 to 
log Li ~ 1, as found previously for NGC 4151. 

Then, using the unabsorbed X-ray luminosity of 
L2-iokeV 4 X 10'*® erg s“® of the highest-fiux epoch as 
a proxy for the ionising luminosity producing O VII and 
the best-fit X-ray spectr al slope of this epoch of F 1.88 
dSteenbrugge et ^l2005f ) to estimate the mean ionising pho¬ 
ton energy, the calculated distance of the coronal line region 
in NGC 5548 from the central ionising source for the above 
gas density is ra 2850 light days 7.8 light years. As 
found previously for NGC 4151, this value places this re¬ 
gion well beyond the hot inner face of the obscuring dusty 
torus, which in NGC 5548 is measur ed by dust reverbera - 
tion campaigns to be ~ 40 — 60 days (iKoshida et al.l[2014l l. 
Furthermore, it is consistent with the fact that a recovery 
of the flux of the coronal lines to the high level seen in the 
2002 epoch is occurring only in 2013 (see Section E21) and 
that only a flux decrease is observed in the 5-yr time period 
spanned by our own data set. 

4.2 The X-ray heated wind scenario and the 

difference between NGC 5548 and NGC 4151 

I Pier fc Void (Il995l l proposed that the coronal line region is 
a layer on the inner part of the dusty torus that becomes 
an efficient coronal line emitter only when evaporated in 
an X-ray heated wind. In this scenario, the wind will have 
undergone adiabatic expansion from its launch location at 
the inner face of the torus until the gas density and distance 
from the ionising source are optimal to give the required 
high ionisation parameter. In the process, the coronal line 
gas will have cooled. In Fig. El we present tentative evidence 
that this scenario applies to NGC 5548, similar to what we 
found previously for NGC 4151. In the bottom panel, we 
compare the temperature of the coronal line gas as measured 
by the line ratio [Fe Vii] A6087/[Fe Vii] A3759 (the higher its 
value, the lower the gas temperature; see Fig. HI) with that of 
the hot dust. We find that the two temperatures behave in 
opposite ways; the temperature of the coronal line gas is high 
when the hot dust temperature is low and vice versa. Only 
the data from April 2002 is an exception to this trend and 
might indicate a reddening event. In the top panel of Fig. 
El we show the continuum flux at rest-frame wavelength of 
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Figure 6. The top and bottom panels show the continuum flux 
at rest-frame wavelength of ~ 1 pm (sampling the accretion disc) 
scaled to the [S III] A9531 line flux and the optical coronal line 
ratio [Fe VIl] A6087/A3759, which is a suitable indicator of the gas 
temperature (the higher its value, the lower the temperature; see 
Fig-HJ, respectively, versus the hot dust blackbody temperature. 


~ 1 /rm, which samples the accretion disc luminosity, versus 
the hot dust temperature. A clear correlation is apparent, 
which indicates that the change in temperature for the hot 
dust is due to direct heating by the central ionising source. 
Therefore, the increased AGN radiation that heats the dusty 
torus also increases th e cooling of the cor onal line gas. 

In the scenario of IPier fc VoitI (Il995h . the dusty clouds 
will be evaporated in an X-ray heated wind more efficiently 
for higher AGN luminosities, which will lead to both an in¬ 
crease in mass outflow rate and a stronger adiabatic expan¬ 
sion. We observe this effect individually for NGC 5548 and 
NGC 4151, but it seems to also apply when comparing the 
two sources with each other. The coronal lines in NGG 5548 
are a factor of ~ 3 — 5 stronger than in NGG 4151 relative 
to the low-ionisation narrow emission lines, and so are both 
the ~ 1 fj,m and ~ 2.1 /rm luminosities (higher by a factor 
of ~ 3 and ~ 2, respectively). This indicates that the X- 
ray heated wind producing the coronal lines in NGG 5548 
is stronger than that in NGG 4151, in which case we also 
expect a stronger adiabatic expansion. The much stronger 
fading of the coronal lines with time observed for NGC 5548 
relative to NGC 4151 is most likely a manifestation of it. 


4.3 Strong coronal line emitters as stellar tidal 
disrnption events 


iKomossa et al.l 20081: 

Gelbord et alJ 20091: Komossa et al. 

2009 

: Wane et al.l 2011 

, 2 OI 2 I: Yane et al.ll2013l: [Rose et al. 


2015fl . A stellar tidal disruption event seems to be the most 


plausible explanation for the strong fading of the coronal 
lines observed in these sources over a time period of sev¬ 
eral years. The tidal disruption of a star by a supermassive 
black hole produces an accretion event, a flare, that sub¬ 
sequently fades on timescales of several months to a year. 
The accretion flare is seen as an X-ray/UV outburst, which 
can produce highly ionised emission lines, including emis¬ 
sion from He II. The X-ray luminosity outburst is expected 
to fade with time as or slig htly shallower than this de¬ 

pending on the typ e of the star llReeslll988l: IPhinnevI[19891: 
iLodato et al.luOOgj ') and subsequently a strong fading of the 
high-ionisation lines is expected. 

Could the strong coronal line variability that we observe 
for NGC 5548 have been caused by a stellar tidal disrup¬ 
tion event? Such an event was raised as a possibility for the 
strong increase (by a factor of > 10) of the He II A4686 emis¬ 
sion line flux observed for this source be tween January and 
May 1984 by IPeterson fc Ferlandl (Il986l ~) and interpreted as 
the addition of a new gas component unrelated to the BELR. 
The recently discovered strong coronal line emitters have 
line ratios [Fe VIl] A6087/[O III] A5007 ~ 0.1 — 1 and have 
[Fe x] emission comparable to or stronger than [Fe VIl] be¬ 
fore the fading of the lines. In NGC 5548, we measure line 
ratios [Fe vii] A6087/[O ill] A5007 ~ 0.03 — 0.07, which puts 
this source in the vicinity of the strong coronal line emitters. 
However, the [Fe x] A6375 emission line is a factor of ~ 5 
weaker than the [Fe VIl] A6087 emission line. As expected for 
stellar tidal disruption events, the high-state of the coronal 
lines in 2002 is accompanied by a high state of the X-ray flux, 
which subsequently fades strongly. The observed X-ray flux 
decrease of a factor of ~ 8 in 3 years roughly obeys the ex¬ 
pected decay law for these events (albeit based on only 
two points). However, the flux of both the X-ray continuum 
and O VII f X-ray coronal line are observed to have recov¬ 
ered by 2013 to their original high state (see Section [3.211 . 
Finally, the He II emission line did not show the continuous 
flux decline observed for most of the coronal lines. Its flux 
varied significantly during the observing period, but fluctu¬ 
ated between increase and decrease from epoch to epoch. 
Given these considerations, it is therefore unlikely that the 
strong coronal line variability that we observe for NGC 5548 
was caused by a stellar tidal disruption event. 


5 SUMMARY AND CONCLUSIONS 

We have presented the second extensive study of the coro¬ 
nal line variability in an AGN. Our data set for the well- 
studied source NGC 5548 is unprecedented in that it in¬ 
cludes five epochs of quasi-simultaneous optical and near- 
IR spectroscopy spanning a period of 5 years and three 
epochs of X-ray spectroscopy overlapping in time with it. 
Our main results are as follows. 


Such a strong variability of the coronal lines, i.e. fading of the 
flux, as we observe for NG C 5548 has been reported so far 
for only one AGN (IC 3599: [&upe et al]|l995l : [Brandt et al.l 
Il995^ and is usually associated with a new class of non-active 
galaxies, the so-called strong coronal line emitters. Several 
of these sources have recently been detected in the SDSS 


(i) Whereas the broad emission lines changed only little 
(by ~ 30 — 40%) and the hot dust emission changed mod¬ 
erately (by ~ 70%), the coronal lines varied strongly (by 
factors of ~ 2 — 4). However, the observed high variability 
amplitude of the coronal lines is mainly due to a flux de¬ 
crease. In fact, none of the coronal lines reacted with a flux 
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increase to the increase in ionising flux observed between 
the May 2004 and January 2006 epochs as the Pa/3 broad 
line and hot dust emission did. 

(ii) We have applied plasma diagnostics to the optical 
[Fe VIl] and X-ray O VII emission lines in order to constrain 
the gas number density, temperature and ionisation param¬ 
eter of the coronal line region. We find that this gas has 
a relatively low density of Ue ^ 10® cm“® and requires a 
relatively high ionisation parameter of log t/ ~ 1, similar to 
our previous results for the source NGC 4151. We estimate 
the distance of the coronal line region in NGC 5548 from 
the central ionising source for the above gas density to be 
rit ~ 7.8 light years. This value puts this region well beyond 
the hot inner face of the obs curing dusty torus (of ~ 2 light 
months; iKoshida et al.l[2014l ~). 

(iii) The relatively high ionisation parameter at large 
distances make it likely that the coronal line region is an in¬ 
dependent entity. One po ssibility is that it is an X-ray heated 
wind as first proposed bv IPier &: VoitI (Il995ll . As in the case 
of NGC 4151, we find support for this scenario in the form 
of a temperature anti-correlation between the coronal line 
gas and hot dust, which indicates that the increased AGN 
radiation that heats the dusty torus appears to increase the 
cooling efficiency of the coronal line gas, most likely due 
to a stronger adiabatic expansion. The strong coronal line 
variability of NGC 5548 can also be explained within this 
picture; relative to the low-ionisation narrow emission lines 
the coronal lines in NGC 5548 are a factor of ~ 3 — 5 stronger 
than those in NGC 4151 indicating a much stronger wind, 
in which case a stronger adiabatic expansion of the gas and 
so fading of the line emission is expected. 

(iv) Since NGC 5548 can be classified as a strong coro¬ 
nal line emitter, we investigate if the strong coronal line vari¬ 
ability could have been caused by a stellar tidal disruption 
event, as is often proposed for this kind of sources. Given 
that the flux of both the X-ray continuum and O VII f X- 
ray coronal line are observed to have recovered by 2013 to 
their original high state and that the He II emission line 
did not show the continuous flux decline observed for most 
of the coronal lines, we consider this scenario an unlikely 
explanation for the strong coronal line variability. 
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